Semiclassical approach to fidelity amplitude 
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The fidelity amplitude is a quantity of paramount importance in echo type experiments. We 
use semiclassical theory to study the average fidelity amplitude for quantum chaotic systems under 
external perturbation. We explain analytically two extreme cases: the random dynamics limit 
-attained approximately by strongly chaotic systems- and the random perturbation limit, which 
shows a Lyapunov decay. Numerical simulations help us bridge the gap between both extreme cases. 

PACS numbers: 05.45.Mt,05.45.Pq,03.65.Sq 



I. INTRODUCTION 

Irreversibility in quantum mechanics cannot be ex- 
plained in the same way as in classical mechanics. Instead 
of measuring how difficult it is to invert the individual 
trajectories, Peres proposed [l[ as a way to understand 
quantum irreversibility to measure the difficulty to in- 
vert the dynamics. To put it in other words, to measure 
how difficult it is to implement a certain Hamiltonian or 
to know the specifics of the dynamics, given that com- 
plete isolation is impossible to attain. As a consequence 
the fidelity -also known as Loschmidt echo (LE) [2[- was 
proposed as a measure of instability of a given Hamilto- 
nian under external perturbations. The LE is commonly 
defined as 



M(t) = \0(t)f 



where 



0(t) = (Vole 



iH e t/h —iHot/h 



(1) 



(2) 



is the fidelity amplitude (FA) and H and H e differ by 
a perturbation which is usually taken to be an additive 
term like eV(q,p), so e is the perturbation strength. The 
importance of the LE lies not only in that it can be used 
to identify quantum chaos but also in that it is a mea- 
surable quantity. 

There is a vast amount of work attempting to char- 
acterize the decay regimes of the LE as universal 0- 
Q. After a short time transient, chaotic systems decay 
exponentially. For small perturbation strength - Fermi 
Golden Rule regime - the decay rate depends quadrati- 
cally on the perturbation strength. For larger perturba- 
tion the decay rate is expected to be perturbation inde- 
pendent, and given by the largest Lyapunov exponent of 
the corresponding classical system. For regular systems 



the behavior is fundamentally different: for small per- 
turbation strength the decay is Gaussian - which can in 
fact lead to faster decay than chaotic systems @,0| ~ an d 
for larger values of the perturbation the decay is power 
law (sec reviews 3, [f|). However, for chaotic systems, 
recent works [cH12j have reported non-universal oscilla- 
tory behavior in the decay rate of the LE as a function of 
the perturbation strength. In the way to elucidating the 
origin of such "anomalous" behavior of the LE one key 
quantity to understand is the FA. The FA is important 
in itself for a number of reasons, the most important be- 
ing that in some of the echo experiments it is in fact the 
quantity that is measured [13l - ll8l ]. 

It is important to remark that in previous works the 
characterization of decay regimes is done on the average 
behavior of the LE. Unlike most papers in the field - 
with the exception of fl9l ] - in this work we concentrate 
on the average FA. The fundamental difference with [l9| 
is that we obtain the decay regimes directly of the aver- 
age FA, not requiring to compute the average fidelity. In 
order to do so we use the semiclassical theory known as 
dephasing representation (DR) (l9rJ2l| . We consider two 
limiting situations and obtain rigorous analytical expres- 
sions for the decay rate in both cases. In the first case we 
suppose the dynamics is completely random. For maps 
on the torus each step is given by drawing two random 
numbers from a uniform distribution. This corresponds 
to the limit of infinite Lyapunov exponent. In this limit 
we obtain an expression for the decay of the FA that is 
valid for all times - up to the saturation point. In fact 
it should be observed for very short times for any map - 
our expression is exact for t = 1. As the Lyapunov expo- 
nent is increased the random dynamics decay regime is 
observed for increasingly longer times. We have already 



presented related results in [22[ . The other limiting case 
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presupposes that the perturbation is completely random. 
That is, after each step of the map the perturbation con- 
tributes with a random phase to each trajectory. It shall 
be seen that this assumption can be justified in the large 
perturbation strength limit because semiclassical calcu- 
lations involve complex exponentials of such strengths. 
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In that case -using the DR and transfer matrix theory- 
we show that the asymptotic decay rate of the FA is con- 
trolled by the largest classical Lyapunov exponent A. In 
other words, we are capable of predicting , analytically, 
the appearance of a Lyapunov regime in the average FA 
rather that the fidelity itself (i.e. the square of the mod- 
ulus of the FA). 

Even though in a realistic scenario one may not see 
clearly the appearance of the limiting behaviours, traces 
of one (or both) of them are usually present. For short 
times, on average the decay rate is that predicted for 
random dynamics. This result is valid for longer times as 
the dynamics becomes effectively more random -i.e. as A 
grows. Then, for large perturbations, there is a crossover 
to an asymptotic regime where the decay rate is A/2, 
independently of the perturbation. The crossover time is 
strongly perturbation dependent. We exhibit numerical 
simulations introducing a perturbation model where the 
amount of randomness can be varied. In this model the 
decay of the FA is well described by a sum of two terms, 
each one corresponding to one of the limiting behaviors. 
The relative weight of the terms depend both on A and 
on the perturbation (amplitude, length scale). 

The paper is organized as follows. Section [TT| gives a 
brief introduction to the dephasing representation (DR). 
This semiclassical approximation is the basis of two an- 
alytically tractable models to be developed in the subse- 
quent sections. The random-dynamics approach, which 
is valid for large Lyapunov exponents, is described in 
Sect. Mil Using the baker-map family as a model of 
chaotic dynamics, and employing the transfer matrix 
method, we argue in Sec. HVI that for large perturbations 
the asymptotic decay is ruled by the Lyapunov exponent. 

Section [V] presents a brief study of mixed regimes, 
where the decay of the FA is bi-exponential. We show 
in some numerical examples that the decay rates are still 
given by the models mentioned above. An empirical ex- 
pression for the crossover time is found. 

We conclude in Sec. I VII with a summary of our results 
and some final remarks. 



II. SEMICLASSICAL DEPHASING 
REPRESENTATION OF THE FIDELITY 
AMPLITUDE 

Recently [L9U2H the dephasing representation was in- 
troduced to give a compact and efficient way to compute 
semiclassically the FA. The derivation of the DR starts by 
replacing the quantum propagators by the semiclassical 
Van Vleck propagator -this is the standard semiclassical 
approach. The first innovation is to use the initial value 
representation [2^, [24| for the Van Vleck propagator. In 
this way a full semiclassical -so-called uniform- expres- 
sion for 0(t) can be obtained. However, this expression 
is still too difficult to be computed and a further approx- 
imation is needed which involves using trajectories of H e 
with slightly different initial conditions but which remain 



close to the trajectories of Hq up to a certain time. The 
validity of this "dephasing trajectories" argument was 
justified using the shadowing theorem One of the 

forms of the FA obtained using the DR looks like 

DB (t) = J dqdpW^{q,p)cxp{~lAS e {q,p,t)/h) (3) 

where W(q,p) is the Wigner function of the initial state 
tp and 

AS £ (q,p,t) = -e f dTV(q(T),p(T)) (4) 
Jo 

is the action difference evaluated along the unperturbed 
classical trajectory. In this way the decay can be at- 
tributed to the dephasing produced by the perturbation 
of the actions -thus the name DR. However numerics 
show that Eq. ([3]) also accounts for decay due to classical 
overlaps (2lj . 

The DR is a useful tool to compute fidelity and to as- 
sess fidelity decay. Actually, when referring to fidelity de- 
cay regimes one refers to an average behavior. Through- 
out this work by average we mean average over initial 
states. The average FA over a number n r of initial states 
is 

I r 

DR {t) = —^2 dqdpW 4 , t (q,p)cxp(-iAS £ (q,p,t)/h), 
Ur i=i J 

(5) 

If in particular, {ipi}™l Q is some complete set, the Wigner 
function of an incoherent sum of such a set is a constant. 
We then get 



DR (t) = - J dqdpcxp(-lAS E (q,p 7 t)/h) 7 (6) 

where V is the volume of phase space. This expression 
was used in (2f| to compute the FA averaged over initial 
states. Henceforth we take V = 1. 

Equation ^ is the starting point of our studies. By 
making some additional assumptions on the dynamics 
and/or the perturbation wc shall be able to derive ana- 
lytical expressions for the decay rate of the average FA 
(next sections) . Being at the basis of our future develop- 
ments, it is necessary to be sure that ([6]) gives indeed an 
accurate description of the FA decay. The DR has suc- 
cessfully passed many tests in the last years 0, [U [2(| , 
especially in quantum chemistry [27l - l30| . There has been 
a similar approach to the dephasing representation in the 
case of linear response functions for electronic spectra 
[3l| . We want, though, to verify its performance in the 
specific systems that shall be used in this paper. 

For the numerical comparisons we have preferred to 
use quantum maps on the torus. These maps possess 
all the essential ingredients of the chaotic dynamics and 
are, at the same time, extremely simple from a numeri- 
cal point view -both at the classical and at the quantum 
levels. The periodic boundary conditions that the torus 



geometry imposes translates in a discretization of both q 
and p upon quantization. The Hilbcrt space is then finite 
dimensional (dimension N) corresponding to an adimcn- 
sional Planck's constant h = l/(2nN) (we assume the 
torus has unit area). Position and momentum basis are 
related by the discrete Fourier transform (DFT) . In such 
a setting the quantum map corresponds to a unitary op- 
erator U |32| . In particular we choose quantum maps 
that can be represented in the split-operator form: 



U 



-i2irNT(p) -i2irNV(q) 



(7) 



The reason of using maps with this structure is twofold: 
(i) Many well known classical maps fall in this category, 
e.g., the standard map, the kicked Harper map, the saw- 
tooth map, as well as some cat maps j32l [33|. (ii) They 
allow a very efficient numerical implementation due to 
the possibility of using the discrete Fourier transform. 
The classical version of U is 



P 



P 



dV(q) 

do 
dT(p) 



(mod 1). 



(8) 



Here we shall consider the perturbed cat maps given by 
34 



P = P-aq + Kf{q) 
q = q — bp + Kh(p) 



(mod 1). 



(9) 



In particular, for the numerics in this section and in 
Sec. Mil we consider the perturbing "forces" 

f(q) = 2tt [cos(27rg) - cos(47rg)] , (10) 
h(p) = 0. (11) 

In Eq. (|9|) we set K <C 1 and a, b such that the resulting 
map is completely chaotic (hyperbolic). In such a case 
the Lyapunov exponent is approximately 



A w log 



1 + ab + y/ab(A + ab) 



(12) 



Time is discrete in maps, so from now on we use the 
integer n to count time steps. We shall use the in- 
dependence of the quantum map U to define the fidelity 
amplitude as 

0(n) = <M(£4) n Wo>. (13) 

Accordingly, the perturbation strength is given by 

e = \K- K'\ . (14) 

In order to calculate the semiclassical DR expression 
for the average FA we define a n s — n q x n p grid of initial 
positions and momenta. Thus we get 



OdhW = ^^2exp[-i2nNAS e (q,p,n)] 



(15) 




6 8 10 12 14 16 18 



FIG. 1. (color online) Average fidelity amplitude as a func- 
tion of time for the perturbed cat map of Eq. (|9]) for a = b = 1 
(A = 0.96; • and x), a = b = 2 (A = 1.76; A and *) 
and a = b = 20 (A = 5.99; ■ and +). Both x, *, and + 
symbols correspond to the semiclassical approximation (DR) 
while •, A and ■ correspond to full quantum evolution. 
The three dashed (grey) lines represent the exponential decay 
exp(— At/2). The solid (red) line is exp(— Tt) with V computed 
semiclassically (see Sec. Illlfl . We set s/H = 6.433 which corre- 
sponds to F = 1.484. The horizontal lines indicate saturation 
values: dashed (black) line is 1/y/ru, with n s = 2 x 10 9 (DR); 



solid (grey) line is 1/yiV x n r with N = 2 (Hilbert space 
dimension) and n r = 8000 (number of initial states). 



where q 



i/riq, p = j/rip (with i = 0, 



j = 0, . . . , rip — 1), and f2 = l/n q n p 
ence is 



1. 



The action diffcr- 



AS s (q,p,n) 



(16) 



with F(q) = J f(q)dq and the points (qi,Pi) are the suc- 
cessive iterates of the classical map of Eq. ([9]) . 

In figure Q] we show the time evolution of the fi- 
delity amplitude, comparing the semiclassical approxi- 
mation (DR) with the full quantum evolution. Three 
perturbed cat maps were considered (with Lyapunov ex- 
ponents A = 0.96, 1.76, 5.99) under the same perturba- 
tion. For the full quantum calculation we used a dimen- 
sion N = 2 18 and averaged over n r = 8000 initial coher- 
ent states. For the DR simulation the integral <j6j) was 
discretized, the total number of initial conditions being 
n s = 2 x 10 9 . 

We see that the agreement between quantum and semi- 
classical calculations is quite good in both cases. This 
makes us confident on the dephasing representation as a 
starting point for the analytical understanding of the av- 
erage FA. In the next sections we shall describe two mod- 
els that explain satisfactorily the main features observed 
in figure [1] for a large Lyapunov exponent the dynamics 
can be considered as being essentially random and the 
decay is single-exponential with a rate that depends only 
on the perturbation (Sect. IIIip . If the Lyapunov expo- 
nent is small and/or the perturbation is large enough the 
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asymptotic decay rate is independent of the perturbation, 
equal to half the Lyapunov exponent, (Sect. ITVf . 

To conclude this section we briefly address the mat- 
ter of saturation values. The average of n s , modulo one, 
complex numbers with uniform random phases goes like 
Therefore for chaotic systems and long enough 

times the DR yields approximately O nT! (t) ~ On 
the other hand, it is well known |35j that the saturation 
value of Eq. (fT5|) is determined by 1/ \/~N, with TV the size 
of the Hilbert space. Therefore in the full quantum sim- 
ulation the saturation value is \/N x n r . Notice that in 
the simulations we chose n s ,n r , N so that the saturation 
values of both quantum and semiclassical calculations ap- 
proximately coincide. 

III. FIDELITY AMPLITUDE FOR RANDOM 
DYNAMICS 




2 4 6 8 10 12 14 
e/h 

FIG. 2. (color online) Decay rate V computed from Eqs. (|19p 
and (|20[) as a function of the rescaled perturbation. The 
points indicate the values chosen in figure [T] (e/h = 6.433) 
and figure [3] (e/h = 3.635). The dashed horizontal lines are 
the corresponding A/2 used in figure 1 (0.96/2, 1.76/2 and 
5.99/2). 



Recently [36| the DR was used to derive analytically 
the decay rate of the FA for perturbations acting on a 
small portion of phase space. The advantage of using 
local perturbations is that we can suppose that the tra- 
jectory becomes uncorrelated between successive hits to 
the perturbed region. It can be argued that the same ef- 
fect can be achieved for non-local perturbations provided 
that A is very large a f22j . 

In this section we study the behavior of Odr(ti) in the 
A — > oo limit by assuming that the dynamics is purely 
random. This evolution is completely stochastic in the 
sense that there is no correlation for the different times of 
the evolution. To compute Oonin), we make a partition 
of the phase space in N c cells. We consider that the 
probability to jump from cell to any other in phase space 
is uniform. Therefore it is straightforward to show that 
the mean FA results 

CW") = 4z E - E exp[-i(AS eJl + ... + AS £ , jn )/h] 



— ^expHASej/Tl) 



(17) 



is the action difference evaluated in the cell 
jk at time k. The continuous limit is approached when 



where AS £Jk 



N c — > oo. So, OoR(n) results 



D R.(n) 



J cxp [— i AS e (q,p)/ ft] dqdp 



(18) 



where AS £ (q,p) is the action difference after one step of 
the map. This exponential decay can be rewritten as 



D R(n) = exp(-rn) . 



(19) 



with 



r = -log J cxp[-iAS 6 (q,p)/h}dqdp . (20) 




FIG. 3. (color online) Same as figure [T] with e/h = 3.635 
corresponding to F = 4.448. 



Equation (|20|) is one of the key results of this work: if the 
dynamics is (uncorrelated) random then the average fi- 
delity amplitude decays exponentially with decay rate T. 
In figure [2] we plot T as a function of the rescaled pertur- 
bation strength e/h for the map of Eq. ([9|). The circles 
mark the perturbation strength values used in Figs.[T]and 
[3] Both in Figs. [T] and [3] the agreement with the random- 
dynamics approach is well observed for initial times. If A 
is large enough (A = 5.99 for ■ and + symbols in Figs.Q] 
and [3]) , the dynamics can be considered as random for 
practical purposes and the decay rate -up to saturation- 
is given by T. In Sec. IIVI we shall see in a specific deter- 
ministic model (the baker family) how a large Lyapunov 
exponent leads to an effective decorrelation of the succes- 
sive points in a trajectory, and therefore the map behaves 
as if it were random. 

We point out that in Figs. Q] and 02 after the inital 
decay following the random dynamics predictions, there 
is a crossover to a different decay rate. In figure [1] the 
second decay rate is A/ 2 while in figure [3] it is not. In 
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Sees. IIVI and [V] we analyze in more detail this crossover 
behavior. 



IV. THE RANDOM-PERTURBATION LIMIT 

In the previous section we showed that a large Lya- 
punov exponent leads to a decay that is well described 
by the random dynamical model. Here we argue that for 
large enough perturbations a Lyapunov decay must be 
expected. 

Key to our analysis is the implementation of the trans- 
fer matrix method to calculate the semiclassical average 
fidelity amplitude. The transfer matrix method is a stan- 
dard tool in Statistical Mechanics for the calculation of 
the partition function of Ising lattices [3714391 ]. Due to 
the analogy between Ising partition functions (for imagi- 
nary temperature) and semiclassical periodic orbit sums, 
the method has also found applications in the Quantum 
Chaos domain |UE3. 

We start by analysing the particularly simple case of 
the baker map and a random perturbation. Then, we 
provide analytical and numerical evidence showing that 
such a Lyapunov decay should also be observed in a larger 
class of chaotic maps submitted to nonrandom (large 
enough) perturbations. 



A. The baker map and the transfer matrix method 

The baker's transformation is an area preserving, 
piecewise-linear map of the unit square defined by [44j 



qi = 2q - (jlo , 
Pi = (Po + Mo)/ 2 , 



(21) 
(22) 



where = [2<?o], the integer part of 2qo. This map 
admits a very useful description in terms of a complete 
symbolic dynamics [45| . A one to one correspondence 
between phase space coordinates and binary sequences, 

(p, q) O . . . fi-2fi-i ■ M0M1M2 • • • , Mi = 0)1) (23) 

can be constructed in such a way that the action of the 
map is conjugated to a shift map. The symbols are as- 
signed as follows: is set to zero (one) when the i— 
th iteration of (p, q) falls to the left (right) of the line 
q = 1/2, i.e. [2q%] = fXi. Reciprocally, given an itinerary 
.../i_2/U-i • M0M1M2 • • •) the related phase point is ob- 
tained through the especially simple binary expansions 



oi+l 



i=0 



2 z+l 



•M0M1 



•M-1M-2 ■ 



(24) 
(25) 



For simplicity, let us initially consider a perturbation 
that depends only on q. Using the binary expansions 



defined above, the DR average fidelity amplitude (|15I16|) 
reads: 



Jo 



ie[F(.poA'iM2---)+-F(-AiiM2/J3---) + --- + -F(-A'n-Ui,iMn + !■■■)] 



(26) 



Here and in the rest of this section, for the sake of a 
lighter notation, we write just e in place of e/K. 

In order to introduce the transfer matrix method it is 
necessary to truncate the binary expansions at a finite 
length L. i.e., 



L-l 



Oi+l ' 



(27) 



Equivalently one may think that we are treating exactly 
a perturbation which is constant over vertical strips of 
equal width 1/2 L (we adopt this point of view). With 
this proviso the integral (|26|) becomes a finite sum: 



a 



2 -(n+L-l) ^ 

fJ,Q,...,[l n + L-2=0,l 



e ie[F(.IJ,o---flL-l)+F(-Hl---HL)+---+F(-Hn.-l—IJ-n+L-2)] ^8) 

The prefactor 2~'" +i ~ 1 ) represents the area (weight, 
probability) of the region in phase space corresponding 
to the finite symbolic string fi , . . . , /i„ + ^_ 2 (the area de- 
pends on the length only, not on the particular code). 
These regions taken together form a partition of phase 
space into disjoint elements. All the trajectories starting 
in a given element have the same action. 

The trick now is to express the sum above as a product 
of n matrices. Define the indices 



ki 



2 L 1 X .flo . ../ii_ 2 , 



k n — 2 
Then we can write 



L-l 



X .fi T 



■ HL+n-2 



feo,...,fc, 



M koM M klM 



(29) 
(30) 
(31) 
(32) 



i,k n 



(33) 



where the nonzero elements of the matrix M are given 
by 



1 



o ieF(.fi —fit-i) 



(34) 



(the convenience of having introduced a factor s/2 will 
become clear soon). It is implicit in this definition that 
ko and ki must satisfy the shift condition, i.e., they must 
share the bits fix . . . [1l-2 [see Eqs. (|29l30p ]. otherwise the 
matrix element is set to zero. The shift condition deter- 
mines that nonzero elements are located on two staircases 
of slope 1/2 ("baker structure"), explicitly given by 



fci = 2k - 2 L Vo + HL-2 



(35) 
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FIG. 4. The transfer matrix M (left) can be split into two 
unitary matrices Ui, U2 (center, right). Nonzero elements are 
shown as black squares (L = 5, schematic). 



(see figure [4} . 

If one defines the unit-norm vector 

|1) = 7 =(1,1,...,1), (36) 
then Eq. (|33[) can be written in the compact way 

0^>) = 2-»/ 2 (l|A/"|l)- (37) 

This is a remarkable result: the decay properties of the 
average fidelity amplitude are embodied in the spectrum 
of a finite matrix. In particular, the asymptotic decay is 
ruled by the largest eigenvalue (in modulus) of M. The 
special staircase structure of the transfer matrix allows 
us to deduce some general properties of its spectrum. 

For a fine discretization, the phase in Eq. (|34|) varies 
by small jumps along the staircases of M, the size of the 
jumps depending both on the characteristic spatial scale 
of F(q) ("correlation length") and on the perturbation 
amplitude e. As the correlation length decreases and/or 
the amplitude grows, neighboring matrix elements be- 
come increasingly decorrelated. This leads us to consider, 
as a reference, an ensemble of random transfer matrices 
having the same structure as M, but whose phases are 
i.i.d. random variables chosen from a uniform distribu- 
tion in [0, 27r]. 

Random or not, because of its particular structure, the 
matrix M can be split as a sum of two unitary matrices 
Ux,U 2 (see figure [p: 

M=-j= (t/i + U 2 ) . (38) 

When U 1 and U2 are random unitary matrices drawn in- 
dependently from the Circular Unitary Ensemble (CUE) 
[32L |4d| . analytical and numerical studies say that the 
spectrum of M lies almost entirely inside the unit circle 
in the complex plane. As the dimension of Ui is increased, 
the largest eigenvalue (in modulus) tends to the unit cir- 
cle. In the limit of infinite dimensional matrices, all the 
spectrum is contained in the unit disk [H, H3| ■ 

We have not attempted an analytical study of the spec- 
tral properties of random matrices with the baker struc- 
ture. However, our numerical calculations show that they 
possess very similar statistical properties to those of the 



normalized sum of two CUE matrices (as far as the spec- 
trum close to the unit circle is concerned). Figure [5] shows 
clearly that most eigenvalues arc distributed inside the 
unit circle, the relative number of "outsiders" decreasing 
as dimension in increased. 




FIG. 5. (Color online) Eigenvalues Zj of random matrices 
having the baker structure. Shown are the combined spectra 
of 64/8 matrices of dimension 64x64/512x512 (left/right). 
Circles of unit radii were drawn for reference. 



So, as the leading eigenvalues of M are close to the unit 
circle, we conclude that for a random perturbation the 
asymptotic decay of O dr (n) is governed by the prefactor 
in (|3"?1) . thus we have 

|0^>)| 2 ~2-". (39) 

This is an exponential decay, with a rate given by the 
Lyapunov exponent of the baker map. This is the other 
key result of the present work: starting from the FA, 
and using the DR and the transfer matrix method we 
obtain that the asymptotic decay of the average FA is 
given by A/2 (consequently the fidelity decays with the 
rate lambda). 

After having analysed the random transfer-matrix 
model, we question to what extent a nonrandom per- 
turbation may lead to a random-matrix behavior. As 
mentioned before, e.g., a smooth but large perturbation 
may cause a randomization of the phases (|34[) . In or- 
der to check this expectation we consider now the baker 
dynamics with the action difference of Eqs. (|16|10p . i.e., 
F(q) = sin(27rg) — sin(4.7rg)/2. The spectrum of the cor- 
responding M{e) is shown in figure [H There we see that 
for, say, e > 10 the leading eigenvalues settle down in 
a small neighborhoo d of the unit circle. Then, for such 
large perturbations |O dr (?i)| 2 decays with a rate which 
is close to the Lyapunov exponent of the classical map. 

Note, however, that in typical numerical simulations 
the decay rate may exhibit some departure from the Lya- 
punov exponent. This may be due to two reasons. First, 
the perturbation amplitude may be not large enough 
for the settling down of the Lyapunov regime. Second, 
even for completely random perturbations, if M is finite- 
dimensional then the leading eigenvalue will be in general 
outside the unit circle. It is true that, as the dimension 
N of M is increased, the leading eigenvalue tends to the 
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Zi 




FIG. 6. (Color online) Transfer matrix eigenvalues (moduli) 
as a function of perturbation amplitude. The perturbation 
F(q) = sin(27rg) — sin(47rg)/2 was discretized in a lattice of 
256 points (L = 7). 



unit circle, but it does so very slowly, possibly like TV" 1 / 3 
[43l ]. So, in practice, some finite deviations from the Lya- 
punov decay rate should not be ruled out. 



B. Extensions 

Here we present some extensions of the application of 
the transfer matrix method to account for more general 
maps and/or perturbations. 

First note that perturbations that depend both on q 
and p can easily be accomodated in the transfer-matrix 
scheme, for, according to (|23|) . it requires only adding 
some symbols for the discretization in p. 

Now consider the i?-adic baker map defined by ([48}) 



Qi = Bq - vq , 
Pi = (Po + vo)/B , 



(40) 
(41) 



where B is an integer > 2, and v = 0,1,..., B — 1. 
This map is completely analogous to the usual (B = 2) 
baker map, but with expansion rate B. At the s ym bolic 
level it corresponds to a full shift on B symbols [45|. Its 
Lyapunov exponent is log B. Thus, by making B large 
enough we shall be able to enter the random-dynamics 
regime (Sect. ITTT|) . 

The transfer matrix M is constructed in a similar way 
as we did for the baker: basically, all we have to do is 
to change from base-2 to base-i? in the standard-baker 
formulas. By making the appropriate substitutions, and 
skipping the details, one arrives at the following expres- 
sion for the average FA 



B (n)^B-^ 2 (l\M^\l) 



(42) 



with |1) a unit-norm vector of length B 1 ^ 1 . 

For a perturbation that only depends on q the non-null 
elements of Mb read 



( M ^) fco , fcl 



1 



JeF(.v ...v L - 1 _) 



(43) 



where Vi are the "decimal" digits of the £>-ary expansion 
of the coordinate q, truncated at depth L. The non-null 
elements of Mb lie along B flights of stairs of slope 1/B 
(riser=l, tread=B). In analogy to the B = 2 case, we 
can split Mb into a sum of unitary matrices: 



M F 



1 B 

Vb f-r 



(44) 



Again, in the case that Ui are large independent CUE 
matrices, it is possible to show that the spectrum of Mb 
is almost completely contained in the unit disk in the 
complex plane |47j. The same behavior is expected for 
the transfer matrix of the baker family, when a nonran- 
dom but large enough perturbation is considered. In such 
a case, one would have the asymptotic decay 



\0 B (n)\~e 



(45) 



with Xb = log-B, the Lyapunov exponent of the classi- 
cal _B-baker map. As an example, we show in figure 
the transfer matrix spectrum for B = 13 and a smooth 
perturbation. Plot are the decay rates 



logS 



- log|zi 



(46) 



where {z{\ stand for the eigenvalues of Mb- 

Focusing on the asympotic decay, we see that for small 
perturbation amplitudes the decay is very well described 
by random-dynamical model, as was expected from the 
discussion of Sect. IIIII Note, however, that the extension 
of the random-dynamic regime is sensitive to the form of 
the perturbation. For large e the leading decay rate fluc- 
tuates around the Lyapunov value (log B)/2, in a similar 
way as we had seen in the case of the standard baker. 

The baker family permits a transparent estimation of 
the crossover point between random-dynamics and Lya- 
punov regimes. For a g-dependent perturbation and a 
B-baker map, the average fidelity amplitude is given by 



O dr (h) = / dq e 
Jo 



ie[F(q )+F(q 1 )+...+F(q n - 1 )] 



(47) 



Here it is implicit that qi , q 2 , . . . are related to q a by the 
dynamics. Alternatively we can take the ^ to be inde- 
pendent and write 



Je[F(q )+F(q 1 )+...F(q n . 1 )] 



DR (n) = J dq dqi ... dq n ~\e %e 
P(q Q ,qi, . . . ,q n -l) , (48) 
where 

P(. . .) = 5( qi -g(qo))S(q 2 -g 2 (qo)) . . . d^-g^ 1 (q )) , 

(49) 

and g(q) = Bq - [Bq]. 

If the joint probability distribution P were factorable, 
then we would recover the formula for random dynam- 
ics. Strictly speaking this factorization is impossible for a 
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FIG. 7. (Color online) Transfer matrix decay rates for the B = 
13 baker as a function of perturbation amplitude. We con- 
sidered two perturbed actions F(q) — sin(27rg) — sin(47r(/)/2 
(top panel) and F(q) — sin(27rg) (bottom panel). These func- 
tions were discretized in a lattice of 2197 points (L = 3). 
Dashed horizontal lines correspond to Lyapunov decay, i.e., 
(log_B)/2. Dotted lines are the predictions of the random- 
dynamics model. 
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<7o <7o 

FIG. 8. The probability distribution is a product of delta 
functions S(qi — g(qo))3(q2 — 9 2 {lo)) ■ ■ ■■ Gridlines define the 
critical square binning. All bins contain the same probability, 
meaning that the histogram is constant, therefore separable. 
The same is true for any bin size larger than 1/B (B = 3 in 
these plots). 



novelty is that some zeros will appear along the staircases 
whenever the associated symbolic sequence is prohibited 
(the information about allowed and prohibited sequences 
is contained in the transition matrix fioT]). The corre- 
spondence between symbols and phase space points will 
be nontrivial, but it could be settled numerically, just by 
launching initial conditions and registering the sequence 
of Markov regions that are visited. 

The main obstacle we see at the moment is that in gen- 
eral the symbolic sequences do not have the same weight. 
It is possible that these weights can be embodied in the 
transfer matrix, or it may be a good approximation to as- 
sume that for long sequences the weights are equal. This 
issue, which exceeds the original scope of the present pa- 
per, is currently being investigated [491 ] . 



V. CROSSOVER REGIME 



deterministic dynamics, however, when considering P in 
the context of the integral above one sees that P admits 
some coarse graining without affecting the integral. This 
is true whenever the scale for the smoothing is smaller 
than the typical scale of variation of the phase exp(ie . . .), 
which is of the order of 1/e. The smaller e, the larger the 
smoothing of P we are allowed to do. Suppose that the 
smoothing consists of substituting P by a histogram of 
bin size a ~ 1/e for each coordinate (see figure [8]). Then, 
by choosing a > 1/B the smoothed distribution becomes 
a constant, i.e., becomes separable (the marginals are 
constant). We conclude that for e < B the distribution 
P is effectively separable and the random-dynamic result 
is applicable. 

Concerning the application of the transfer matrix 
method to more general maps with a finite symbolic dy- 
namics, e.g., cat maps, it is likely that the scheme can 
also be adapted to this case. The general structure of 
the transfer matrix will be the same as that of the corre- 
sponding baker with the same number of symbols. The 



In the previous sections we provided analytical descrip- 
tions of the FA in two extreme cases: the random dynam- 
ics limit and the random perturbation limit. The first 
model provides a satisfactory description for initial times 
in general and for all times - before saturation - in the 
large A limit (ideally A — > oo). The second explains the 
appearance of the Lyapunov regime for long times and 
large enough perturbations. 

In many cases the time behavior of the FA is not 
pure, but shows a crossover from a random-dynamics de- 
cay (short times) to a random-perturbation decay (long 
times). We know that for short times the influence of the 
dynamics is not significant on average, then the trajecto- 
ries in the DR formula (|15[) behave like random walks 
with uncorrelated jumps. Therefore, for short times 
the decay of the average fidelity amplitude is given by 
Eq. -it is exact for t — 1. In addition, for a gen- 
eral smooth perturbation with large amplitude, we have 
shown that there is a randomization of phases yielding an 
effectively random perturbation. As a consequence, the 
well known Lyapunov decay emerges in the large e limit. 
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There remains a need for understanding the crossover (in 
the time domain) between random-dynamics and Lya- 
punov decays. 

In order to investigate the crossover we resorted to nu- 
merical simulations using the DR and the cat map ([9]) 
with a random perturbation. We divided the torus into 
N c x N c square cells, assigning a random constant value 
vi € [-0.5,0.5] to each cell, with i e [1,N%]. Next we 
computed DR (n) by summing over n s initial conditions 
and then averaged |0 DR (n)| over N IC realizations of the 
perturbation. 

In figure [9] we show the results of our simulations for 
various values of e and A. Using Eqs. (fT9|) and (|20f we 
computed the decay rate T for the initial times (see inset). 
We can already see on the top panel that < |0 DR (n)| > 
has three well defined regimes: two exponential decays 
followed by the expected saturation due to finite number 
of initial conditions. Then, up to the saturation point, 
the decay is well described by 



< |0 DR (n)| >~exp(-Tn) + ,4exp(-An/2). (50) 

We found empirically that, in our particular model, the 
prefactor A is given by 



A = X/N c . 



(51) 



So, when A — > oo the dynamics is random, the second 
term vanishes, and the decay rate is T for all times before 
saturation, as predicted in Sec. Mil When A is finite but 
the perturbation is completely random (N c — > oo) the 
first term in Eq. (|50|) dominates the initial decay. The 
reason for this is that deterministic motion together with 
random perturbations looks, on short time scales, like 
random dynamics (as far as the accumulation of phases 
is concerned). 

In figure [9] (top and middle panels) the dependence 
of A with A is clearly observed. The dependence with 
N c is esta blished in the bottom panel, where we plot 
< |0 DR (n)| > for three different orders of magnitude of 
the randomness parameter N c . We remark that, as ex- 
pected, the crossover point depends both on the dynam- 
ics (through A) and on the perturbation (through N c ). 
We can clearly see this in the middle panel of figure [9] 
where we show < |0 DR (n)| > for different values of the 
rescalcd perturbation strength e/h [and for A = 0.96, cor- 
responding to a = b = 1 in Eq. and in the bottom 
panel of the same figure where we plot < |0 DR (n)| > for 
one value of e/h and N c = 10, 100 and 1000. We can 
again see that Eq. (|50|) describes quite well the decay up 
to the saturation value. 

We thus verify that the crossover regime can be de- 
scribed by a natural bi-exponential law combining the 
random-dynamic and Lyapunov decays. 



VI. CONCLUSIONS 

We have analyzed the average fidelity amplitude, 
a quantity that characterizes the instability and irre- 
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FIG. 9. (color online) Top: < |0 DR (rc)| > for e/h = 6.2589 (+ 
symbol in the inset), and different values a = b = 1, 2, 3 and 
4 in Eq. © giving (o) A = 0.96, (A) A = 1.76 , (o) A = 2.389 
and (V) A = 2.887 respectively. Middle: < |0 DR (n)| > for 
a = b = 1, A = 0.96 and different perturbation strengths: (+) 
e/h = 6.2589; (x) e/h = 6.794; (*) e/h = 8.235 ; (□) e/h = 
43.318. Inset: T(t) computed from N TC ~ 100 realizations 
of the perturbation with N c = 100. The points indicate the 
value of e/h for the curve with corresponding symbol on the 
bottom panel. Bottom: < |0 DR (ra)| > for e/h = 6.2589 (+ 
symbol in the inset), for three different values of N c : (•) 
iV c = 10, (A) N c = 100 and (?) iV c = 1000 (a = 6 = 1, 
A = 0.96). The dashed blue line (in all panels) corresponds 
to the bi-exponential in Eq. (|50[) . The horizontal gray lines 
correspond to the expected saturation values. 



versibility of perturbed quantum evolution. We have 
taken as a starting point a well established scmiclassi- 
cal theory -the dephasing representation Our first key 
result is obtained in the limit of random dynamics. We 
show that in that case the decay of the average FA is 
exponential for all times - up to the saturation - and we 
give an analytic expression for the corresponding decay 
rate. Maps on the torus with a very large Lyapunov ex- 
ponent are expected to mimic random dynamics rather 
precisely. Thus the decay rate T obtained in equation 
PU)) is to be observed for increasingly larger periods of 
time as A is increased. Numerical results shown in figures 
1 and 3 confirm this assertion. 

Our second key result is that using the DR and a novel 
approach in this context - the transfer matrix method - 
we were able to predict that the asymptotic decay of the 
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average fidelity amplitude should be perturbation inde- 
pendent and given by A/ 2. We are even able to explain 
possible deviations from this behavior. We point out that 
the Lyapunov, as well as other regimes were obtained in 
[l9j analyzing statistics of action differences in an ex- 
pression which is obtained from the square of the aver- 
age FA (equation (O). However in that article the short 
time behavior due to random dynamics-like behavior is 
not observed. In addition, one major advantage of our 
approach is that we do not need to sum over pairs of tra- 
jectories - a difficulty that arises in the squaring needed 
to perform average fidelity studies. 

As consequence, we have both analytical and numerical 
evidence to conclude that the average FA for an arbitrary 
chaotic map, should decay initially with a decay rate 
given by the random dynamics approximation, and for 
long times there is the Lyapunov regime. We cannot yet 
predict in a generic case the behavior or the time where 
the crossover takes place. Nevertheless in an attempt to 
understand this problem we have studied a model of ran- 
dom perturbations which helped us in determining that 
for random perturbations and intermediate times the de- 
cay is also given by the Lyapunov exponent. 

In [22| it is shown that the decay of Loschmidt echo has 
three regimes as a function of the perturbation strength. 
The well known -quadratic- Fermi golden rule regime 



is valid for small perturbations. Then there is a regime 
where the decay rate is given by T (computed for ran- 
dom dynamics) and after that there is a crossover to a 
perturbation independent regime given by the Lyapunov 
exponent A. We believe that this paper is an important 
step towards understanding this complex behavior of the 
LE. We remark that the average FA is tightly related to 
the LE but has some fundamental differences. In fact, 
if we suppose that n s is the number of initial conditions 
then the average LE can be expressed as a sum of the 
square of the average of the FA (times n s ) minus a sum 
of "nondiagonal" terms consisting of products of 0(t) for 
different initial conditions. Here we have taken a signif- 
icant step forward by understanding one of the parts of 
the LE: the average FA, which is in itself an important 
quantity in many experiments. In order to fully com- 
prehend the behavior of the LE further work is needed 
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